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Abstract 

In this paper, the crucial phenomenon of the expansion of the uni- 
verse has been discussed. For this purpose, we study the vacuum 
solutions of Bianchi types / and V spacetimes in the framework of 
f{R) gravity. In particular, we find two exact solutions in each case 
by using the variation law of Hubble parameter. These solutions corre- 
spond to two models of the universe. The first solution gives a singular 
model while the second solution provides a non-singular model. The 
physical behavior of these models is discussed. Moreover, the function 
of the Ricci scalar is evaluated for both the models in each case. 
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1 Introduction 

The accelerating expansion of the universe has attracted much attention in 
recent years. Why the universe is expanding at an increasing rate and spread- 
ing its contents over greater dimensions of space? It is the most interesting 
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question in the modern day cosmology. There is another issue of dark energy 
and dark matter which has been discussed widely. Einstein gave the concept 
of dark energy by introducing the cosmological constant. Later, he remarked 
that the introduction of the cosmological term was the biggest blunder of his 
life. 

Carmeli pQ predicted for the first time in 1996 that the universe expansion 
was accelerating. Later, the supernova experiments [2] became known and 
it was found that major part of the universe was filled with dark matter 
and dark energy. The results from these experiments provided evidences 
for the accelerating expansion of the universe. The cosmological constant is 
considered as one of the candidate responsible for the dark energy. Modified 
theories of gravity provide an alternative approach to study the universe. 

There are some useful aspects [3] of modified theories of gravity. Mod- 
ified gravity gives an easy unification of early time inflation and late time 
acceleration. It provides a natural gravitational alternative to dark energy. 
The explanation of cosmic acceleration is obtained just by introducing the 
term 1/R which is essential at small curvatures. The modified gravity also 
describes the transition phase of the universe from deceleration to accelera- 
tion. It can be used for the explanation of hierarchy problem in high energy 
physics. The f{R) theory of gravity is one of the modified theories which 
is considered most suitable due to cosmologically important /(-R) models. 
These models consist of higher order curvature invariants as functions of 
the Ricci scalar. Viable f{R) gravity models [1] have been proposed which 
show the unification of early-time inflation and late-time acceleration. The 
problem of dark matter can also be addressed by using viable f{R) gravity 
models. 

Singularity has been an important issue in General Relativity (GR). The 
occurrence of spacetime singularity is a general feature of any cosmological 
model under some reasonable conditions. It may be possible to avoid these 
undesired singularities in the context of modified theories. Kanti et al. [5] 
showed that there did not exist any cosmological singularity by considering 
higher order curvature terms. Nojiri and Odintsov [6j proposed some re- 
alistic singularity free models in modified f{R) gravity. Bamba et al. [7] 
proved that non-minimal gravitational coupling can remove the finite-time 
future singularity in modified gravity. Thus there is a strong reason to study 
solutions of the field equations in modified theories to address the problems 
like dark energy and singularity. 

Weyl [S] and Eddington [9J studied f{R) actions in 1919 and 1922 respec- 
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tively. Buchdahl [TU] explored these actions in the context of non-singular 
oscillating cosmologies. Cognola et al. [Hj investigated f{R) gravity at one- 
loop level in de-Sitter universe. It was found that one-loop effective action 
can be useful for the study of constant curvature black hole nucleation rate. 
Spherical symmetry is the most closest approach to the nature because one 
can compare the results from the solar system observations. Thus the most 
commonly explored exact solutions in f{R) gravity are the spherically sym- 
metric solutions. Multamaki and Vilja [12] studied spherically symmetric 
vacuum solutions for the first time in this theory. The same authors |T3] also 
investigated the perfect fluid solutions and showed that pressure and density 
did not uniquely determine /(-R). Capozziello et al. [H] explored spherically 
symmetric solutions of f{R) theories of gravity via the Noether symmetry 
approach. Hollenstein and Lobo [15] analyzed the exact solutions of static 
spherically symmetric spacetimes in f{R) gravity coupled to non-linear elec- 
trodynamics. 

Cylindrical symmetry is next to spherical symmetry which may be used to 
study the exact solutions of the field equations in f{R) gravity. Azadi et al. 
[16j studied cylindrically symmetric vacuum solutions in metric f{R) theory 
of gravity. Momeni [17] extended this work to the general cylindrically sym- 
metric solution. Recently, we have explored static plane symmetric vacuum 
solutions [H] in f{R) theory of gravity. The field equations are solved using 
the assumption of constant scalar curvature which may be zero or non-zero. 
However, no attempt has been made so far for solutions with non-constant 
scalar curvature. 

The accelerating expansion of the universe can be studied using Bianchi 
types / and V spacetimes which are the generalization of FRW spacetimes. 
Due to the spatially homogeneous and isotropic nature, many authors |19j - 
[21] have studied these spacetimes in different contexts. Berman [25] intro- 
duced a different method to solve the field equations by using the variation 
law of Bubble's parameter. The main feature of the variation law is that it 
gives the constant value of deceleration parameter. Using this law, Singh et 
al. PB] explored perfect fluid solutions of Bianchi type V spacetime. Re- 
cently, Kumar and Singh [27] studied solutions of the field equations in the 
presence of perfect fluid using Bianchi type / spacetime in GR. The same au- 
thors [28] investigated perfect fluid solutions using Bianchi type / spacetime 
in scalar-tensor theory. 

In this paper, we focuss our attention to explore the vacuum solutions of 
Bianchi types / and V spacetimes in f{R) theories of gravity using metric 
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approach. The paper is organized as follows: In section 2, we give a brief 
introduction about the field equations in the context of f{R) gravity. Sections 
3 and 4 are used to find exact vacuum solutions and the singularity analysis 
of these solutions. In the last section, we summarize and conclude the results. 

2 f{R) Gravity Formalism 

The f{R) theory of gravity is the generalization of GR. The action for this 
theory is given by 

Here f{R) is a general function of the Ricci scalar and is the matter 
Lagrangian. It is noted that this action is obtained just by replacing R 
by f{R) in the standard Einstein- Hilbert action. The corresponding field 
equations are found by varying the action with respect to the metric g^i, 

F{R)R^, - ]^f{R)g^u - V^V,F{R) + g^,,UF{R) = kT^,, (2) 

where 

F{R) = df{R)/dR, □ = VV^, (3) 

is the covariant derivative and T^,^ is the standard matter energy-momentum 
tensor derived from the Lagrangian L^. These are the fourth order partial 

differential equations in the metric tensor. The fourth order is due to the 
last two terms on the left hand side of the equation. If we take f{R) = R, 
these equations reduce to the field equations of GR. 
Now contracting the field equations, it follows that 

F{R)R - 2f{R) + 3nF{R) - kT (4) 

and in vacuum, we have 

F{R)R-2f{R) +3nF{R) ^0. (5) 

This gives an important relationship between f{R) and F{R) which will be 
used to simphfy the field equations and to evaluate f{R). 
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3 Exact Bianchi Type / Solutions 



Here we shall find exact solutions of Bianchi I spacetime in f{R) gravity. For 
the sake of simplicity, we take the vacuum field equations. 

3.1 Bianchi Type / Spacetime 

The line element of the Bianchi type / spacetime is given by 

ds"" = dt" - A\t)dx^ - B\t)dy^ - C\t)dz\ (6) 

where A, B and C are cosmic scale factors. The corresponding Ricci scalar 
becomes 

A B C AB BC CA^ 

^ = -2[i + 5 + c + I5 + 5C + cm]' 

where dot represents derivative with respect to t. 
We define the average scale factor a as 

a = \/ABC (8) 

and the volume scale factor as 

V = a^ = ABC. (9) 
We also define the generalized mean Hubble parameter H in the form 

H = ^{H^ + H^ + H^), (10) 

where Hi = ^, H2 = ^, H3 = ^ are the directional Hubble parameters in 
the directions of x, y and z axis respectively. Using Eqs.(l8|)- (fT0|) . we obtain 

H=\y- = l-{H, + H, + H,) = ^. (11) 
6 V 6 a 

It follows from Eq.dS]) that 



Inserting this value of f{R) in the vacuum field equations, we have 

F{R)R^, - V^V,F{R) _ F{R)R - nF{R) 
9t^u 4 



(13) 



Since the metric (E]) depends only on t, one can view Eq. (fT5]) as the set of 
differential equations for F{t), A, B and C . It follows from Eq. fll3l) that the 
combination 

^ F{R)R,, - V,V,F{R) ^^^^ 

is independent of the index fi and hence — Ay = for all fi and u. 
Consequently, Aq — Ai = gives 

B C Ab CA Af F ^ 
-B-C^AB^CA^AF-F = '- ^''^ 
Also, Aq — A2 = and Aq ~ A3 = yield respectively 

0, (16) 

: 0. (17) 

Thus we get three non-linear differential equations with four unknowns namely 
A, B, C and F. The solution of these equations can be found by following 
the approach of Saha [29] . 

3.2 Solution of the Field Equations 

Subtracting Eq.^ from Eq.([l5]), Eq.^ from Eq.^ and Eq.^ from 
Eq. ffT^ . we get respectively 

AO IP A h 
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These equations imply that 

B 

C , ^ f dt 

A , ^ f dt 
^ = rf3exp[c3 



f^.,exp[c,/-^], (21) 
§ = ^2exp[c2 I ^], (22) 
^ c/3exp[c3 [ 4^], (23) 
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where Ci, C2, C3 and di, d2, are constants of integration which satisfy the 
relation 

ci + C2 + C3 = 0, rfirfarfs = 1. (24) 
Using Eqs. (pT]) - (!23|) . we can write the metric functions exphcitly as 

/dt 
(25) 

f dt 

B = ap2exp[q2 / (26) 

/dt 
(27) 



where 
and 



Pi = {dr'd2'')K p2 = id,d2~')K ps = {d,d2^)-^ (28) 

2ci + C2 Ci - C2 Ci + 2C2 

?i = ^ , ?2 = — - — , qs = • (29) 

Notice that pi, p2, Ps and qi, q2, qs also satisfy the relation 

P1P2P3 = 1, qi + q2 + q3 = 0. (30) 

Now we use the power law assumption to solve the integral part in the 
above equations. The power law relation between scale factor and scalar field 
has already been used by Johri and Desikan [30j in the context of Robertson- 
Walker Brans-Dicke models. However, in a recent paper [21], Kotub Uddin 
et al. have established a result in the context of f{R) gravity which shows 
that 

where m is an arbitrary constant. Thus using power law relation between F 
and a, we have 

F = ka"", (31) 

where k is the constant of proportionality and m is any integer. The decel- 
eration parameter q in cosmology is the measure of the cosmic acceleration 
of the universe expansion and is defined as 

da , , 

, = -51- (32) 
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It is mentioned here that the negative sign and the name "deceleration pa- 
rameter" are historical Initially, q was supposed to be positive but recent ob- 
servations from the supernova experiments suggest that it is negative. Thus 
the behavior of the universe models depend upon the sign of q. The positive 
deceleration parameter corresponds to a decelerating model while the nega- 
tive value provides inflation. We also use a well-known relation [32] between 
the average Hubble parameter H and average scale factor a given as 

H = la-"", (33) 

where / > and n > 0. This is an important relation because it gives the 
constant value of the deceleration parameter. 
Using Eqs.([II]) and ([33]), it follows that 

a = la^-"" (34) 
and consequently the deceleration parameter turns out to be 

q = n-l (35) 
which is obviously a constant. Integrating Eq. fl34|) . it follows that 

a={nlt + ki)^, riy^O (36) 

and 

a = k2exp{lt), n = 0, (37) 

where ki and ^2 are constants of integration. Thus we obtain two values 
of the average scale factor that correspond to two different models of the 
universe. 

3.3 Singularity Analysis 

The Riemann tensor is useful to determine whether a singularity is essential 
or coordinate. If the curvature becomes infinite at a certain point, then the 
singularity is essential. By constructing scalars from the Riemann tensor, 
it can be checked whether they become infinite somewhere or not. It is 
obvious that infinitely many scalars can be constructed from the Riemann 
tensor. However, symmetry considerations can be used to show that there 
are only a finite number of independent scalars. All others can be expressed 



8 



in terms of these. In a four-dimensional Riemann spacetime, there are only 
14 independent curvature invariants. Some of these are 



D D _afe p p p pafe p p ryabcd p pafe pcd 

ni — n — g nab-, -K2 — -Kab-K , -K3 — Uabcd-ti , -(14 — -K^^/t^j. 

Here we give analysis for the first invariant commonly known as the Ricci 
scalar. For Bianchi type / spacetime, it is given by Eq.([7]). For a special case 
when m = —2, it follows from Eg. (13 II) that 



F = ka-'^. (38) 
After some manipulations, we can write 

^3P(aa + d^) - (gig 2 + 52^3 + 93^1)1 

which shows that singularity occurs at a = 0. 



^1 = -^i ^27^ J y^^) 



3.4 Model of the Universe when n ^ 0. 

Now we discuss model of the universe when n 7^ 0, i.e., a = {nit + fci)". For 
this model, F becomes 

F = k{nlt + ki)~^ . (40) 

Using this value of F in Eqs. (!25l) - (p71) . the metric coefficients A, B and C 
turn out to be 

.qi{nlt + ki)~^ . 



A = pi{nlt + ki) " exp[- 



kl{n - 1) 



p . , , ^i Mnlt + ki) n 



kl{n-l) 





(41) 


n ^ 1 


(42) 


n ^ 1. 


(43) 



C = Mnlt + k,)-^eM 
The directional Hubble parameters Hi {i = 1, 2, 3) take the form 

H = ^ + (44) 

' nlt + ki k{nlt + ki)^ 

The mean generalized Hubble parameter becomes 

H = , ^ , (45) 
nit + ki ^ ^ 
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while the volume scale factor turns out to be 

V={nH + ki)n. (46) 
Moreover, the function of Ricci scalar, f{R), can be found by using Eq. f|T2|) 

f(R) = ^(nlt + ki)^R + 3kf{n - l){nlt + k^)^^ . (47) 

It follows from Eq.([39]) that 

R^R, = -2\?,l\2-n) {nit + k,r'- ^^^ + ^^3 + ^3^1) ^^^^ + k,p^^^ {AS) 

which clearly indicates that f{R) cannot be explicitly written in terms of R. 
However, by inserting this value of i?, f{R) can be written as a function of 
t, which is true as R depends upon t. For a special case when ^^ = |, f{R) 
turns out to be 



f{R) = ?R 



3fc/2 -9/2 ± y^81/4 + tal±^3±£3Mi? ^ 



(49) 



2 ' 2R ' ' 

This gives f{R) only as a function of R. 

3.5 Model of the Universe when n = 0. 

The average scale factor for this model of the universe is a = k2exp{lt) and 
hence F takes the form 

k 

F = —exp{-2lt). (50) 

Inserting this value of F in Eqs.f l25p -( l271) . the metric coefficients A, B and 
C become 

A = p,fc2exp(/t)exp[- ^^"^P^~^^^ ], (51) 

klk2 

B = p,k2exp{lt)exp[-^-^^^^^], (52) 

KIK2 

n 1 fu\ r ^3exp(-/t) 

C = p3fc2exp(/t)exp[ — J. (53) 
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The directional Hubble parameters Hi and the mean generalized Hubble 
parameter will become 

Hi = l + ^ exp(-/t), H = 1. (54) 
kk2 

The volume scale factor turns out to be 

V = k2^ exp{3lt) (55) 

while f{R) takes the form 

f{R) = J^exp{-2lt){R-QP). (56) 
For this model, R becomes 



,2 (?i?2 + 92^3 + qsQi] 



R = -2[Ql'- ' "t^^n . (57) 

^ A;2A;2^exp(2/t) ^ ^ ' 

Here we can get the general function f{R) in terms of R 

f{R) = [R^ + 6fR - 721^] (58) 

2(gig2 + g2g3 + g3gi)^ ^ ^ ^ 

which corresponds to the general function f{R) |33j . 

/(i?) = ^a„i?", (59) 
where tt, may take the values from negative or positive. 

4 Exact Bianchi Type V Solutions 

Here we shall find exact solutions of the Bianchi type V spacetime in f{R) 
gravity for the vacuum field equations. 

4.1 Bianchi Type V Spacetime 

The metric for the Bianchi type V spacetime is 

ds^ = de - A\t)dx^ - e^'^'=[B\t)dy^ + C^{t)dz\ (60) 
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where A, B and C are cosmic scale factors and m is an arbitrary constant. 
The corresponding Ricci scalar is 



ABC 3m2 AB BC CA, , , 

With the help of Eq.([Tl]), we can write - Ai = 0, - = and 
— ^3 = respectively as 

B C 2m2 AB CA AF F 

+ T7^ + 77T + T^-T^ = 0, (62) 



B C A^ AB CA AF F 
A C 2m2 AB BC BF F 



+ in5 + l^ + irB--B = ^^ (63) 



A C A^ AB BC BF F 
A B 2m2 BC CA CF F 



A B A^ ^ BC^ CA^ CF F °' ^^^^ 
The 01-component can be written by using Eq.(l2]) in the following form 

We solve these equations using the same procedure as for the Bianchi type 
/ solutions. 

4.2 Solution of the Field Equations 



Here we get same equations Eqs. (118 1) -( 1201) as obtained previously with the 
difference of the constraint equations (using Eq. (!65!) ) 

pi = 1, p2= P3'^ = P, qi = 0, q2 = -qs = Q- (66) 
Consequently, the metric functions become 

dt , 



A = a, B = aPexplQ J C = aP-^ exp[-Q J 



i^F^ 



(67) 



Moreover, Ri turns out to be 



^1 = -2 [ ^2 ] (68) 



which also yields a singularity at a = 0. 



12 



4.3 Model of the Universe when n 7^ 0. 



Using the value of F from Eq. (HO|) in Eg. (1671) . the metric coefficients A, B 
and C become 

A = (nlt + h)^, (69) 
B = nnlt + h)^exp[ ^^''^^^^^'l^ ], n^l (70) 



C = p-\nlt + k,t^eM ^^^!{n_^l\ " I (71) 
The directional Hubble parameters ifi, H2 and -ffs take the following form 

H2 = / , + r, (73) 

nlt + ki k{nlt + ki)n 

I Q 
H2 = -r, , r- 74 

nlt + ki k{nlt + ki)n 

Notice that the mean generalized Hubble parameter if, the volume scale fac- 
tor V and /(-R) turn out to be the same as for the Bianchi type I spacetime. 

4.4 Model of the Universe when n = 0. 

For this model, the metric coefficients A, B and C turn out to be 

A = k2exp{lt), (75) 

B = Pfc,exp(/t)exp[- ^"^P^~^^^ ], (76) 

klk2 

C = P-%exp{lt)eM^^^^^]- (77) 
The directional Hubble parameters Hi, H2 and take the form 

H, = l, i/, = ; + «f^, = , - . (78) 

kk2 kk2 

Here we also have the same mean generalized Hubble parameter H, volume 
scale factor V and f{R) as given in Eqs. flM]) - flSB]) . 



n-l 
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5 Summary and Conclusion 



The main purpose of this paper is to discuss the well-known phenomenon 
of the universe expansion in the context of f{R) gravity. For this purpose, 
we have investigated exact solutions of the Bianchi types / and V space- 
times using the vacuum field equations. We have found two exact solutions 
for both spacetimes by using the variation law of Hubble parameter. This 
yields constant value of deceleration parameter. These solutions correspond 
to two models of the universe. The first solution gives a singular model with 
power law expansion and positive deceleration parameter while the second 
solution provides a non-singular model with exponential expansion and neg- 
ative deceleration parameter. It is mentioned here that the solutions for both 
spacetimes correspond to perfect fluid solutions [261 122j in GR. We have also 
evaluated function of the Ricci scalar, f{R), for both models in each case. 
In particular, the general function f{R) includes squared power of the Ricci 
scalar for the non-singular model. The physical behavior of these models is 
given below. 

First we discuss singular model of the universe. This model corresponds 
to n 7^ with average scale factor a = [nit + fci)™ . It has a point singularity 
at t = = — ^. The physical parameters Hi, H2, H-^ and H are all 
infinite at this point but the volume scale factor vanishes here. The function 
of the Ricci scalar, f{R), is also infinite while the metric functions A, B 
and C vanish at this point of singularity. Thus we can conclude from these 
observations that the model starts its expansion with zero volume aX t = ts 
and it continues to expand for < n < 1. 

The non-singular model of the universe corresponds to n = with average 
scale factor a = /c2exp(/t). It is non-singular because exponential function 
is never zero and hence there does not exist any physical singularity for this 
model. The physical parameters Hi, H2, H^ are all finite for all finite values 
of t. The mean generalized Hubble parameter H is constant while f{R) is also 
finite here. The metric functions A, B and C do not vanish for this model. 
The volume scale factor increases exponentially with time which indicates 
that the universe starts its expansion with zero volume from infinite past. 

It would be worthwhile to study other Bianchi type spacetimes especially 
by removing the vacuum condition. It would also be interesting to explore 
Bianchi types / and V solutions with perfect fluid. These are under progress. 
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